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We measure the ground and excited states for B mesons in the static limit using maximally variance reduced 
estimators for light quark propagators. Because of the large number of propagators we are able to measure 
accurately also orbitally excited P, D and F states. We also present some results for A;,. 



1. INTRODUCTION 

Usually one calculates the light quark propa- 
gators needed for hadronic observables in lattice 
QCD by exact inversion from one source. It is 
necessary to iterate almost to machine precision 
to avoid introducing bias. This is clearly waste- 
ful, as the variation from one gauge to another 
can be large. The situation is even more severe 
for static quarks, when one gets only one number 
per inversion. 

Stochastic estimators for an inverse of a posi- 
tive definite matrix M are easily obtained from 



Ml 



(1) 



Unfortunately the Wilson-Dirac fermion matrix 
Q is positive definite only for extremely small val- 
ues of hopping parameter K . This problem can 
be solved by taking M = Q^Q, and calculating 
the propagator from 



G ij = Qji = ((Qik4>k)*fa) 



(2) 



This stochastic estimate can be used to cal- 
culate any hadronic observable, but because one 
is usually interested in exponential decay of cor- 
relators at large T values, the variance coming 
from stochastic inversion will kill the signal rather 
quickly. Therefore one would prefer to use im- 
proved operators; this can be achieved using vari- 
ance reduction. 



2. MAXIMAL 
REDUCTION 



VARIANCE 



One suggestion is to use local multi-hit vari- 
ance reduction for <fi fields jD ; this is very easy to 



implement and provides significant error reduc- 
tion. However, as the action is quadratic in <fi one 
should be able to implement a method that takes 
into account not only the nearest neighbours of 
4>, but all fields inside some given region P. 

If the scalar fields inside region P and on the 
boundary of P are called 4> and s respectively, one 
can obtain maximally variance reduced estimate 
for fax) by 



D<j> 
1 



exp{--<p*M jk <j) k + cj)*M jk s k + h.c.) (3) 

where we have also distinguished the elements of 
M connecting fa to only fields inside P (M) and 
those connecting fa to the boundary (M). The 
integral (0) is Gaussian and one obtains easily 



v t = ~M tj L M ]k s k . 



(4) 



We will call the maximally variance reduced 
stochastic estimator for fa. 

To form propagators, we need a product of two 
4> fields. This means that it is not sufficient to 
improve only one field within P, but one must 
choose two disjoint regions P and R and solve for 
two variance reduced fields v and w respectively. 
Then one can apply (g) to obtain 



Gij = {(QikV k )*Wj) 
with no bias. 



3. B MESONS AND A fc 



(5) 



Consider mesons containing one infinitely 
heavy quark Q and one light quark q. Such a 
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system is obtained in leading order heavy quark 
effective theory HQET § for the B meson. This 
system is particularly suitable for stochastic in- 
version as usually one is only able to obtain one 
number per inversion. 

The expectation value one has to calculate con- 
tains simply one light propagator and the prod- 
uct of gauge links U. This is easily obtained from 
equation (||). In addition one can use smearing 
(or fuzzing) to create different sources; in par- 
ticular it is possible to construct operators corre- 
sponding to orbitally excited states by using tech- 
niques from ||]. 

It is also possible to construct observables with 
more than one light quark. Ah in leading order 
HQET contains one infinitely heavy quark and 
two light quarks. This is also easily obtained from 
stochastically improved fields v, w but a little care 
is needed. One way to grasp the subtlety is to 
imagine that there are two quarks with different 
flavours. Then one has to split the sum over sam- 
ples a = 1 . . . N into subsets for each flavour. If 
these subset are independent, one obtains prop- 
agators of each flavour with no bias. In practice 
if the stochastic estimators <\> a are independent 
of a, one can calculate the required propagators 
from 

G % G i'f = £(Q«^Qw«&)X< ( 6 ) 

The Ab correlator is then easily constructed by 
multiplying two light quark propagators from 
equation (Q) by the gauge links corresponding to 
a heavy quark propagator. 

4. RESULTS 

We have performed simulations with tadpole 
improved quenched clover action (csw = 1-57), 
with f} = 5.7. The lattice size we use is rather 
small, 12 3 x 24. We work with two different hop- 
ping parameters Ki = 0.14077 and K 2 = 0.13843. 
The same values of hopping parameters are used 
also in (J] , and we use those results also to set our 
scale. The values K\ and K% correspond to the 
strange quark and twice the strange light quark 
masses respectively. 

We use 20 gauge configurations, and in each 
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Figure 1. Effective mass plot for 5- wave with 
hopping parameter K\ = 0.14077 corresponding 
to roughly strange quark mass. 



of them we calculate 24 stochastic samples for 
each hopping parameter. The stochastic update 
consisted of 5 overrelaxation steps followed by one 
heat bath. Measurements were taken after every 
25 combined OR/HB sweeps and the system was 
thcrmalised for 50 heat bath sweeps before the 
first measurement. 

The choice of the regions P and R is rather 
arbitrary; we chose to divide the lattice by time 
planes at T = 1 and T = 9. We used two different 
fuzzing levels for all operators. 

A typical effective mass plot can be seen in Fig- 
ure |l|, where we have plotted the effective mass of 
the L = (S) state together with factorising fit. 

Our full results can be seem from Figure |^, 
where we have plotted the mass splittings be- 
tween excited states and ground state in units of 
ro ■ Also included is the mass splitting of A& and 
the ground state of B. Note that the ordering of 
P-waves we observed before || has now changed. 

These results should be compared to results ob- 
tained using usual inversion techniques. In Fig- 
ure b| we have plotted results several other groups 

— |S[] have obtained using much more computing 
resources. Our results are clearly consistent with 
theirs. However, we have much smaller error- 
bars and are able to obtain reliably several ex- 
cited states, which is not generally true for the 
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Figure 2. The spectrum of B meson containing 
one light (with mass proportional to pseudoscalar 
mass m 2 p) and one infinitely heavy quark. Also 
included are our results for A D . 



earlier work in the static limit. One should also 
note that our results are compatible with results 
obtained from non-relativistic QCD |J. 

Because we have not yet estimated systematic 
errors arising from finite lattice size and spacing, 
we chose not to extrapolate our results to the chi- 
ral limit. 

5. CONCLUSIONS 

We have shown that stochastic estimators for 
light quark propagators are extremely useful in 
heavy quark effective theory, where the usual 
method of calculating light quark propagators is 
extremely wasteful. We have obtained reliable 
signals for the first time for excited states in the 
static limit, and show that even states with L = 3 
(F-waves) can be obtained. 

It is obvious that our method can be used suc- 
cessfully not only for spectroscopy but also for 
matrix elements. The ability to obtain reliable in- 
formation on excited states is valuable when one 
is trying to extract the heavy-light decay constant 
1b- 



Figure 3. Comparison with earlier results. FNAL 
group U has only estimated the excited S-wave 
with no error. UKQCD results for A& are from 
, Wuppertal results for PO from || . The exper- 
imental values 0Jll[] are for real B meson and 
A{,, so one would expect to see some difference 
from leading order HQET. 
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